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Abstract 

We study the dynamics of a second order phase transition in a situation that mimics a 
sudden quench to a temperature below the critical temperature in a model with dynamical 
symmetry breaking. In particular we show that the domains of correlated values of the 
condensate grow as ^/t and that this result seems to be largely model independent. 
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I. Introduction 

Most of the studies in finite temperature field theory deal with problems in thermal 
equilibrium. There are systems though, where conditions change so fast that to take the 
equilibrium formalism, even as a first approximation, would be inappropriate. One of 
the situations where non equilibrium phenomena are interesting is the dynamics of phase 
transitions. They have been considered by a number of authors, mostly in connection 
with problems in cosmology. More recently the interest in phase transitions ocurring 
in relativistic heavy ion collisions has grown due to the possibility of producing in the 
laboratory the deconfined phase of strong interactions. Clear signals that the quark gluon 
plasma was indeed produced and the study of its its properties is made dificult by the fact 
that in the laboratory this phase is produced in a small region and is very short lived. Non 
equilibrium phenomena are quite normal ocurrence in this kind of experiment and this 
provides another motivation for considering the dynamical aspects of phase transitions in 
relativistic field theories. 

One very interesting proposal on how to decide whether the quark gluon plasma was 
created or not during a heavy ion collision was suggested many years ago and since then has 
been studied by many authors [1] . It is based on the fact that besides the deconfining phase 
transition, we expect in QCD another phase transition at about the same temperature, 
namely, the chiral symetry restoration. The picture goes as follows. The value of the 
order parameter -^-i/; becomes zero at high temperatures reached in the plasma formed by 
the colliding nuclei and chiral symmetry is restored. With the expansion of the plasma 
though, this temperature decreases and eventually chiral symmetry is broken again. The 
direction of symmetry breaking in isospin space, however, does not necessarily have to be 
the same as it originally was before the collision took place. Therefore, if the value of 
the order parameter is, say, in the tt^ direction, after hadronization this plasma will be 
converted into neutral pions only, and not charged ones. This could lead to spectacular 
events with large number of pions correlated in isospin space coming out of one collision 
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and signaling clearly that chiral symmetry restoration was achieved inside the plasma. 
Events of this kind might have already been seen in cosmic rays (the so called Centauro 
events [2] ) as it has been sugested [3] . The dificulty with this picture is that the direction 
in which chiral symmetry is broken does not to have to be the same all over the plasma. 
More likely there will be domains in which the order parameter is the same, but different 
domains will have diferent symmetry breaking directions. This domain structure is a 
very general phenomenon in condensed matter, cosmology, etc. If the typical domains are 
small and only a few pions can be made with energy contained in each of them, then the 
correlation in isospin space can hardly be distinguished from that of a random production. 
In equilibrium there are indications that this is indeed the case. Recently Rajagopal and 
Wilczek [3] have suggested that if the rapid cooling of the plasma is modelled as a sudden 
quench below the critical temperature the size of correlated domains can be much larger 
than the size in equilibrium, growing with time as some power of time. This conclusion has 
been supported by analogies with condensed matter systems obeying classical dissipative 
stochastic first order equations. However, as has beens stressed by Boyanovsky et al. [4] 
the situation in a relativistic field theory, where the order parameter obeys a quantized 
conservative second order equation does not necessarily have to be the same. The same 
problem of determining the rate of growth of domain structures out of equilibrium is also of 
importance in cosmology, since the size of the domains determine the number of different 
topological defects produced during the phase transition. In this paper, we report on 
an attempt to study this phenomenon systematically within the context of a model that 
exhibits dynamical symmetry breaking. 

The method we will use to analyze this problem is the Closed Time Path formalism 
invented in the sixties by Schwinger and Keldysh [5] . Since chiral symmetry breaking is not 
well understood in QCD even at zero temperature, we will take the Nambu-Jona-Lasinio 
[6] model as an effective theory having the essential physics of chiral symmetry break- 
ing present in QCD. In this model the critical temperature is a function of the coupling 
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constant g and, instead of coupling our system to some heat bath to generate the quench- 
ing, we will mimic the drop in temperature by suddenly changing g, and consequently 
the critical temperature. Ideally one should study such a phenomenon starting from an 
initial temperature T (and coupling constant g such that T > Tc{g) ) and then suddenly 
increasing g to g' at some instant to such that T < Tc{g'). The system, in this case, will 
find itself in the wrong phase and the formation of domains in the broken phase will take 
place. Our aim is to calculate the value of < 'ip{t,'r)'ip{t,'r)'ip{t,0)'ip{t,0) > which gives a 
measure of the distance scale over which isospin is correlated. However, such a calculation 
is quite complicated and to bring out the qualitative features we choose an alternate but 
parallel model where the interaction (g) is turned on suddenly at t = such that the initial 
temperature T < Tdg). We find that the domain size grows with time (for large times) as 
y/i in a model independent manner. 

II. Closed Time Path formalism 

In this section we briefiy review the Closed Time Path formalism. All essential details 
can be found in the references [5,7,8]. Suppose we are interested in the expectation value 
of some observable A at instant t in a given ensemble. It is given, in general by 

<A>{t)^ Trp(t)A = Trp(0)C/(0, t)AU(t, 0), (1) 

where p is the density matrix that describes the (mixed) state of the system and U is the 
evolution operator (in the Schroedinger picture) corresponding to a hamiltonian H that 
is , in general, time dependent. The density matrix does not necessarily have to commute 
with the hamiltonian, in which case it describes a non equilibrium state. For t < we 
assume that the system is in equilibrium so that at t = 0, p can be written as 

P(0) - ^v^, (2) 
for some initial hamiltonian Hi. For positive times we take H as the hamiltonian which 
governs the dynamics of the system. We can write now 

Tre-f^H, TrU{T-il3,Ty ^' 
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where T is any time T < 0, since U{T — if3,T) involves only the hamiltonian Hi. We 
can now go back to equation (1). After substituting equation (3), inserting the identity 
operators 1 = C/(0, T)C/(T, 0), 1 = U{t, T')U{T', t) and rearranging the product inside the 
trace we obtain 

TrUjT - T)U{T, T')U{T', t)AU{t, T) 
TrU{T -il3,T)U{T,T')U{T',T) ' 

where T' is some time larger than t. Equation (4) can be pictured as describing the 
evolution of the system from T to T' (with an insertion of the operator A sX t), and then 
backwards in time to T and finally along the imaginary time axis to T — ij3. We can 
now take the limits T — > — oo and T' — > oo. This suggests the definition of a generating 
functional 

Zc[Jc] = Tr Uj{-oo - i/3, -oo)Uj{-oo, oo)Uj{oo, -oo), (5) 

where the subcript c is to remind us that the quantity is defined in the contour in the 
complex t-plane described above. Uj stands for the evolution operator under the infiuence 
of the external sources Jc- Clearly, if H equals and is time independent, and if the 
source Jc is the same along the entire time contour, Zc is the partition function of the 
system and this formalism reduces to the well known imaginary time formalism. We are, 
however, dealing with a nonequilibrium phenomenon where H{t) is different along the 
time contour and, consequently we choose the external sources to be different along the 
two branches of the contour to allow us to obtain Green's functions by taking derivatives 
with respect to the sources. Zc has a path integral representation 

Zc[Jc] = [ D4> e'i'c(^+-^<^)'^'^ (6) 

J (anti)periodic 

with (f) denoting a generic field of the theory. As usual we can separate the quadratic part 
of the lagrangian and expand the interaction part in a power series to obtain a perturbative 
expansion. It turns out that the free propagators vanish if they connect a point in the 
segment (— oo, — oo — to a point on the real axis. Thus the contribution of this segment 
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effectively decouples and contributes only to an overall normalization [8]. Therefore, we 
are effectively left with an integral on the real axis (both ways). It is convenient to rewrite 
this as a normal integral over the real line along one branch, but the propagator then 
acquires a 2 x 2 matrix structure 



Aa6 = 

where 



A++ A+- 
A_+ A_. 



(7) 



A++(a; 


-y) 


=< T{<t>{x)<t>{y)) >p 


A__(a; 


-y) 


=< T*{cf>{x)cf>iy)) >p 


A+_{x 


-y) 


=< (t>{y)(l){x) >f3 


A_+ix 


-y) 


=< (j){x)(j){y) >f3 . 



(8) 



Here T and T* stand for time ordering and antitime ordering respectively. The matrix 
structure and the different boundary conditions satisfied by the functions Aab can be 
understood remembering that in the part of the contour that runs in the negative direction 
(we call it " — " as opposed to the other one that we call " +" ) the time ordering gets reversed, 
and that any point in the " +" branch is earlier than any point in the " — " branch. Besides 
the matrix structure of the propagator the other difference in relation to the usual Feynman 
rules at T = is that we have the interaction term in both branches of the contour and 
consequently there are two kinds of vertices, one connecting fields living on the "+" branch 
and the other (with opposite sign) connecting fields defined on the " — " branch. 

III. The Model 

The model we would like to consider is the Nambu-Jona-Lasinio model in four dimen- 
sions, with two flavors and N colors. 

C = i^a^d^iJa + -^[(i^ai^af - (V^aTsT^a)'] , (9) 

where a = 1,...,N, t*'s are the generators of flavour SU{2), and the flavor indices are 
suppressed for simplicity. This model is nonrenormalizable in four dimensions. However, 
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we treat this, for our discussions, as an effective low energy theory for energy scales below 
the cut off A. The theory has the chiral symmetry 

(10) 

As usual we can introduce auxiliary fields cr and tt* to write the Lagrangian (9) in an 
equivalent form 

As is well known, the chiral symmetry (10) of the theory is dynamically broken if > 27r^ 
and restored at = const. (1 — ^)^A. In our discussion we will take the coupling 
constant g to be time dependent with the simple form 

g{t) = e{t)g, (12) 

and the initial temperature T smaller than the critical temperature 

277^ 1 

Tc = const. (1 - -^)2A. (13) 

Thus, according to the model we have chosen, for t < we have a gas of massless free 
fermions in equilibrium at a temperature T and the state of the system is invariant under 
the chiral transformation of equation (10). At t = we turn on the interaction such that the 
system suddenly finds itself at a temperature much below the critical temperature given 
in (13). This guarantees that the effect of the sudden change will make the chiral symmetric 
initial state of the system unstable, almost in parallel to the case of an expanding plasma 
as discussed in the introduction. An important point to note, however, is that, since both 
the dynamics and the initial state are chirally symmetric, a will have a vanishing vacuum 
expectation value in the ensemble. However, in each element of the ensemble, different 
regions of the space will have different directions of symmetry breaking and, although the 



expectation value of a (averaged over the ensemble) remains zero, correlations of the a 
(and tt) field will grow. 

The complete u ( or tt ) propagator is obtained from 

Dai i^^y) = ~\^,y) + ^ab{x, y) , (14) 

where a, 6 = ±, the zeroth order propagator and S is the self-energy. The discussion of 
domain growth, however, is best carried out in terms of the physical propagators (namely, 
the retarded, advanced and correlated ones) and, therefore, we transform from the ± basis 
to the physical basis through the transformation [5,7], 

D' = VDV~^, S' = VT.V~^, (15) 



with 



Thus, 



\/2 VI 1 



^'=i V. ^n'^V (18) 



with 

Dr = D++ - D+_ 
Da = -D++ + D_+ 
= D++ + 

= + S+_ 

= — — E — 1_ 
Ec = + S 

The functions -Dit(A) are the usual retarded (advanced) functions defined as 

Dr{x, y) = ^(x" - yO) < [0(x), 0(1/)] >p 
Da{x, y) = 9{y^ - x^) < [0(a;), 0(y)] >p, 



(19) 



(20) 



and which gives the fermionic correlation fuunction we are interested in is 

D,{x,y)=<[<^{x),^{y)]+>p . (21) 

In case of Fermi fields we would exchange commutators by anticommutators and vice versa. 

The difference in sign in the definition of cind D^j^j^^ is due to the fact that S_| 

and E |_ include a "— " sign coming from a "— " type vertex. Using equation (14) we find 

the relation between the propagator and the self energy to be 

and 

Dc = DrEcDa + n-- (24) 

Note that in a non equilibrium situation these functions will depend on each time argument 
separately, and not only on their difference. In other words, in a non equilibrium situation 
there is no time translation invariance. Thus, we can not diagonalise them by taking a 
Fourier transform and the order of the functions in (22, 23, 24) is important. Consequently, 
unlike the equilibrium case, it is in general difficult to invert expressions such as 1 + EijD^. 

IV. The Calculation 

The functions I)it(A) carry information about the dynamics of the model and at zeroth 
order are independent of the initial state (they are expectation values of c-numbers which 
vanish at equal times by microcausality). It is the function D^. that contains the infor- 
mation about the state of the system, like particle number distributions, temperature, etc. 
For instance, for free fermions we have (We use the metric r/^i/ = (H )) 



Sc(x) = j ^e-'' '(ki'j^ + m) 2m(2n(ko) - me - m^), 



± ieko 



(25) 



where 



is the Fermi distribution function. We want to calculate D,, for the a (and tt) propagators 
at equal times. At tree level these are not real (on shell) propagating particles, consequently 
they do not thermalise. Another way of saying this is that they are too heavy (their masses 
are of the order of the cut-off) and they are Boltzmann suppressed at any temperature 
much smaller than the cut-off. We have then 

iArA2' 

nO _ (27) 
Dl = 0. 

We will calculate the self-energy in the leading order in so that the main contribution 
comes from the fermionic one loop bubble diagram. But we note first that since g{t) — 
0{t)g, ^R(A){xjy) Q'lid T,c{x,y) vanish if either < or < 0. Consequently, they have 
the generic form 

S(x,|/) = ^(xO)E(x,|/)^(|/0). (28) 

where the quantity with an overbar is calculated with a constant g much the same way as 
in an equilibrium calculation. As a result, the inversion of quantities such as (l-hSijD^)"^ 
becomes trivial for positive times. For example, (1 -|- T,jiD^)~^ satisfies the equation 

{1 + ^rD%)-\1 + ErD%) = 1. (29) 

Written out explicitly, this has the form 

J d'z (1 + 1:rDI)-\x, z)[5\z - y) + ^e{z')i:a{z - y)e{y^)] = 5\x - y). (30) 

By definition Sij vanishes for z^ < y^ and consequently, if we restrict to > 0, the two 
step functions in the second term can be ignored and the inversion can be done much the 
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same way as for a constant coupling, equilibrium case. The knowledge of (1 + ^ for 

positive times is enough to calculate the retarded propagator for positive times (a;°, y° > 0) 

DR{x,y)=[dz D%{x,z){l + J:RD%)-\z,y) 

J (31) 

= Dr{x - y). 

The same argument goes through for the advanced function as well but this is not always 
the case. For the Feynman Greens function, for example, the analogue of equation (30) 
is a complicated integral equation. This comes about because S4.+ describes propagation 
in both, forward and backward time direction so, even after restricting the final points 
x^ and j/° to be positive, 2;°, the intermediate time coordinate, can be negative. Since 
the dynamics for negative times is differentfrom that at positive times the inversion is not 
straightforward. Since our non equilibrium retarded and advanced functions are the same 
as the ones in equilibrium (always for x^,y^ > only) we have (It will become clear shortly 
that this is all we need for our discussion.) 

^^'-■^' = /(S-~ -.^A^U.W '''' 

From equations (24) and (27) then, we obtain for x^,y^ > 

Dc{x,y)= fd^zd^z' DR{x,z)J:a{z,z')DA{z',y) 

r (34) 

= j d^zd^z' DR{x-z)e{z'')i:^{z-z')e{z''')DA{z' -y). 

Note that to evaluate Dc we need Dr(^a) only for positive times as stated earlier. To find 
TiR, T,A and Sc we first calculate E06, a,b — ± using the diagrammatic methods and then 
use the definitions in (19). S^b is simply obtained from the Feynman rules discussed above 
to be 

Safe = -ab g^N Tr j iSab{k)iSUP + k), (35) 
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where Sab ^ire the finite temperature, real time propagators at the initial temperature T 
S++{p) = {-f^'pi, + m) ( ^ ^2—^ + 27Tin{po)6ip^ - m^)] 

\P ~ ITT' + / 

S—{P) = (t'^Pp + m) (^- ^2 -m^-ie + 27rin(po)5(p' - m^)^ ^^^^ 
S+-{p) = {i^p^ + m)2TTi{n{po) - e{-po))5{p^ - m^) 



2 2\ 

m 



S-+{p) = (7''p^ + m)27ri(n(po) - 0{po))5{p 
and "a6" is a sign factor coming from the fact that the "-" vertices carry a negative sign. 

Note that in our theory m = and using these definitions the calculation of flab is 

straighforward. We only give the results here: 

E++=i4:g N / ——3 { 2 [pQUk-p.k) 

S{{po + Ukf - a;p+fe)(^sflfn(po) + ri{ujk) - n{ujk)n{ujp+k))] 
+ Pq^ -po] 



{n{ijjk) - l)(n(po + i^k) - 0{-Po - <^k)) 
+ (-Po - P-k)S{{po - ujkf - i^l+k) 

n{u)k){n{po -ujk) - 9{-pQ+uj))] 
d^k 1 

^2^ — [(^^o^fe - P-k)5{{.Po + ^kf - ^l+k) 



(37) 



n{uJk){n{po - ujk) - 0{-po + Uk)) 
+ (-Po - P-k)S{{po - uJkf - '^l+k) 

{n{iOk) - l)(n(po + <^fc) - 0{-pQ - Uk))]. 
We can now construct Si?, T^a and Sc- They are given, after some rearrangement and 
with the help of the change of variables A; — > — A; — p in various places (these are finite 
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integrals since the momenta are cut off at A), by 

T.R{p) = iAg N / — — ) — -— —+PQ,e^-pQ,-e\ 

J {2'ny 2uk {po -i^k+ ley - uj^^^ (38) 

and 

/d'^k 1 1 
— [(^((Po + uJk? - i^l^j,){-sgn{po) + n{uk) ^^^^ 

- n{u}k)n{po + ujk)) +Po Po]- 

Consequently we can also obtain DfK^j^^ in a straightforward manner. We will need to know 
the position of the pole(s) ep(e*) of Dii(^a){poiP) in order to evaluate D^. We determine 
the position of the pole in the following way. At the pole of Da we have 

ziVA2 + EH(po,p) = 0. (40) 

If we set p = in the above expression, we obtain 

A. fl _ + |lpSln(Zl^) - '-^ r * 1 ^ = 0. (41) 

V 27r2y Stt^"^" pI Jo e^'' + 1 pi - Ak^ ^ ' 

A simple graphical analysis shows that, for > 2n^ (the symmetry breaking case) and 
temperatures smaller than Tg, there is always a solution of equation (41) withpg negative. 
Let us call this solution — M^. To find the next term in the series expansion we substitute 
Pq — > — + bp^ and keep terms up to second order in p. After some tedious algebra we 
obtain 

il-b)ls + {l-b)h = 0, (42) 



where 



and 



/3 = AM' £ dk tgHPk/2) (^^^4lP)3 (43) 
h = IQM^ £ dk tgK/5V2)p^4W" ^^^^ 
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The solution to equation (42) can be easily seen to lie within the range 

l<i<l. (45) 

The position of the pole up to this order is given by Pq = = — + with b 
determined as above. This way we have 

6^-~±z(M-^). (46) 

The wrong sign of the mass (— M^) is a signal of instability of the initial state. One could 
go on to higher orders but this will suffice for our discussion. 

Using our previous results we have now, with x — {t,r) and y = {t,0) 

D,{x,y)= [ d^zd^z'e{zy{z'^) [ 5t /f -ik.ix-z)-ip.(z-z')-n.iz'-y) 
J ^ ^ ^ ^ J (27r)4 (27r)4 (27r)4 

X DR{k)i:c{p)DA{q) ^^^^ 

J 27r 27r h - Po - po - qo + ie 

X Dii{ko,k)T,c{Po,k)DA{qo,k). 

Note that it is the limited range of integrations on z^, z'^ that produces the denominators 
shown above instead of the usual 8 functions that assure conservation of energy. As the 
system is unstable, T)r grows exponentially with time so that one should be careful in 
defining its Fourier transform. It is well known [9] that if ~ e"* for large t, the Fourier 
transform Dji{u!) will be analytic only above the line Im co = a and the inverse Fourier 
transform is given by an integral along a line just above Im co = a. The opposite happens 
for the advanced function, namely, the region of analyticity is Im lu < a. With this in 
mind we can perform the integrations over ko and go to obtain 

^e-^^'^'Dniko, k)- = - zJ«(6fc, k)- ie-^^o^Dfl(po, k) 



/ 



^e^^ 'DA{qo, k) ^ =idA{el, k)- + ie^^^'DA{po. k), 

27r go - Po - «e e - Po 
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where (iR(^) are the residues of -D^(^) at the pole efe(efc)- As is clear from the earlier dis- 
cussion, efc is complex leading to both exponentially growing as well as damped behaviour. 
For large times, it is the exponentially growing solution which will dominate. Furthermore, 
the terms not involving ek do not grow with time and are the ones that would be there 
had g been kept constant and the range of integration in z^,z'^ was (—00, 00) (namely, for 
the case of equilibrium) . We will omit this term from now on since the exponential growth 
gives the dominant contribution. Putting this back in (47) we have 

Dc{t,r-)= [ ^^e^^^-|J.i(6,,^)pe-^(--^^)^^-^S,(po,^) 



(27r) 



where 

f{k ,T)=g N ———— [ I ; ——^i--n{up) + n{up)n{ujp+k)) 

+ ^^7^ — ^p+k) y _ ^p_^^) _ ri{u!p) + n{u!p)n{u}p+k)) 

iLOp — LOp+k — Cfcl 2 

+ — '^^P+k) M 1^ 1 ^^^^^ _ _ n{ujp) + n{ujp)n(ujp+k)) 

((Wp + Wp+fc)2 -^2) 1 / ^ . / X / 

\U!p + iOpj^k — ^k\ Z 

(50) 

A little analysis shows that the function / vanishes for /c^ — > so, for small k"^ we can 
write {k = \k\) 

f{k',T)^T^g{k/T), (51) 

with g{x) as a; ^ as a power, e.g. g{x) ~ x*^, a > 0. We evaluate D^. in equation 
(49) for large times using the saddle point approximation. The saddle point is given by 

k. = (62) 

Therefore, A; — > as t — > 00 and we are justified in disregarding terms of order 0(p^) in 
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the expansion of in equation (46) 



This gives for large times 



Dc{t,r) 



C{T) (^)"e^V 




C{T) (9%^^e 



(53) 
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where C{T) is a function independent of (t, r) and L(t) is the typical size of the domains 
given by 



This is the same scaling as obtained by Boyanovsky et al. [4] for the scalar 0^ theory with 
a mass square that suddenly changes sign. The size of the domains, of course, depends on 
the details of the theory since it is determined by h and M but the power law dependence 
on t in equation (54) seems to depend only on the fact that the first non constant term 
in the expansion of ep for small p is quadratic in p. We see from (53) that the correlation 
grows indefinetely with time but this is an artifact of our approximation. This growth will 
actually stop when inside a domain the value of cr is that of the equilibrium and the phase 
transition will be over. This however would involve a rearrangement of the chiralities of 
the fermions that can not be achieved in leading order in since in this case the fermions 
do not scatter. 

V. Conclusion 

We have shown how to use the Closed Time Path formalism to describe the dynamics 
of a phase transition. In particular we analysed the dynamics of domain growth following 
quenching in a simple model exhibiting chiral symmetry breaking. We found that the 
domains with correlated isospin grow as and we have argued that this result is not 
dependent on details of our model. 
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